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PACS 13 . 85 .Hd - Inelastic scattering: many-particle final states 
PACS 13 . 85 . Lg - Total cross sections 

Abstract - The TOTEM experiment has measured the charged particle pseudorapidity density 
dN c h/d»7 in pp collisions at y^s = 7 TeV for 5.3 < \n\ < 6.4 in events with at least one charged 
particle with transverse momentum above 40 MeV/c in this pseudorapidity range. This extends 
the analogous measurement performed by the other LHC experiments to the previously unexplored 
forward n region. The measurement refers to more than 99% of non-diffractive processes and to 
single and double diffractive processes with diffractive masses above ~3.4 GeV/c 2 , corresponding 
to about 95% of the total inelastic cross-section. The dN c h/dr; has been found to decrease with |?j|, 
from 3.84 ± O.Ol(stat) ± 0.37(syst) at = 5.375 to 2.38 ± O.Ol(stat) ± 0.21(syst) at \n\ = 6.375. 
Several MC generators have been compared to data; none of them has been found to fully describe 
the measurement. 
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Introduction. — The pseudorapidity density of 
charged particles produced in high energy proton-proton 
(pp) collisions reflects the strong interaction dynamics 
that is only partly described by perturbative QCD. Non- 
pcrturbative models and parametrisations are used in the 
Monte Carlo (MC) event generators to describe the hadro- 
nisation of the partonic final states and to model diffrac- 
tive processes [TJ[2]- In the forward region, where periph- 
eral diffractive processes are important, the uncertainties 
are pronounced. A better understanding of these effects is 
also important for the interpretation of high energy show- 
ers recorded by cosmic ray experiments [3H5]- A direct 
measurement of the forward particle densities is, there- 
fore, extremely valuable in constraining the theoretical 
models for particle production in pp interactions. The 
measurement of the charged particle pseudorapidity den- 
sity (dN c h/d7y) in the range 5.3 < \r/\ < 6.4 is presented 
here. This quantity is defined as the mean number of 
charged particles per single pp collision and unit of pseu- 
dorapidity 77, where 77 = — ln[tan(#/2)], and 9 is the polar 
angle of the direction of the particle with respect to the 
counterclockwise beam direction. 

Experimental apparatus. — TOTEM is a dedicated 
experiment to measure the total cross section, elastic scat- 
tering and diffractive processes at the LHC [51 [7] • The 
experimental apparatus [5], composed of three subdetec- 
tors (Roman Pots (RP), Tl and T2 telescopes), is placed 
symmetrically on both sides of Interaction Point (IP) 5, 
shared with the CMS experiment. All three subdetec- 
tors have trigger capability. The Roman Pot stations, 
equipped with silicon detectors and placed at 147 and 220 
m from the IP, detect elastically and diffractively scat- 
tered protons with a small scattering angle down to a few 
pi ad. The Tl and T2 telescopes, placed at about 8 and 
14 m from the IP respectively, detect charged particles 
produced in the polar angular range of a few mrad to 
-100 mrad. The Tl telescope (3.1 < \rf\ < 4.7) consists 
of Cathode Strip Chambers, while the T2 telescope (5.3 
< |?7| < 6.5) is made of triple-GEM (Gas Electron Multi- 
pliers) chambers [9j- The present analysis is based on mea- 
surements with the T2 detector that consists of 2 quarters 
with 10 semicircular chambers each, on both sides of the 
IP. Each chamber provides two-dimensional information of 
the track position in an azimuthal coverage of 192 with a 
small overlap region along the vertical axis between cham- 
bers of two neighbouring quarters (TDJ . Every chamber has 
a double layered read-out board containing two columns 
of 256 concentric strips (400 pm pitch, 80 pm width) for 
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the measurement of the radial coordinate and a matrix of 
1560 pads, each one covering A77 x A<fr ~ 0.06 x 0.018 rad, 
for the measurement of the azimuthal coordinate and for 
triggering. Radial and azimuthal coordinate resolutions 
are about 110 pm and 1 °, respectively jTTj . The total ma- 
terial of 10 chambers amounts only to — 0.05 Xo [12] . The 
read-out of all TOTEM detectors is based on the "VFAT" 
front-end ASIC, which provides as output a digital signal 
and trigger [P3] . 

Detector simulation. - The TOTEM software [13] . 
based on the CMS framework [TS], embeds the necessary 
interfaces to the GEANT4 [TFJ simulation toolkit, the de- 
scription of the material placed between the IP and T2 
and the offline software used for the event reconstruction. 
The description of the T2 detector is implemented within 
this framework. The GEM signal digitisation has been 
parametrised using a dedicated model that, after proper 
tuning, reproduces well the measured cluster (a group of 
neighbouring strip or pad channels) size and reconstruc- 
tion efficiency as a function of the ionisation energy re- 
leased in the gas by the incident particle, diffusion coeffi- 
cient of the fill gas, chamber gain and the VFAT thresholds 
[THH7]. A special effort was devoted to understand and 
quantify secondary particles produced by the interaction 
of particles with the material in front of and around T2 
and then seen in the detector. The simulation of the for- 
ward region, properly tuned with the data, showed that a 
large number of secondary particles is produced in the vac- 
uum chamber walls in front of the detector, in the beam 
pipe conical section at \rj\ = 5.53 and at the lower edge of 
the CMS Hadron Forward (HF) calorimeter [18]. 

Track reconstruction and alignment. The track 
reconstruction is based on a Kalman Filter-like algorithm 
|17| that is simplified due to the small amount of material 
traversed by the particle crossing the 10 GEM planes and 
to the low local magnetic field in the T2 region. The parti- 
cle trajectory can, therefore, be successfully reconstructed 
using a straight line fit. The reconstructed tracks have at 
least 4 hits (pad clusters with or without an overlapping 
strip cluster), of which at least three have a pad/strip clus- 
ter overlap. A x 2 -probability greater than 1% is required 
for the straight line fit. 

Using a coordinate system with the origin located at the 
nominal collision point, the X axis pointing towards the 
centre of the LHC ring, the Y axis pointing upward (per- 
pendicular to the LHC plane), and the Z axis along the 
counterclockwise beam direction, two geometrical track 
parameters are defined: Zo and ZImpact. Zo is the Z value 
at the position of the minimum approach of the track to 
the Z axis; the ZImpact parameter (sec fig. [1]) is the Z 
coordinate of the intersection point between the track and 
a plane ("^2") containing the Z axis and orthogonal to the 
plane defined by the Z axis and the track entry point in 
T2 ("tti"). Due to the short lever arm of the T2 detector 
(~ 40 cm) compared to the distance to the IP, the ZImpact 
and Zo resolution is of the order of 1 m. 
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Fig. 1: Definition of the track ZImpact parameter. The de- 
scription is discussed in the text. 

The relative position of the detector planes within a T2 
quarter (internal alignment) and the overall alignment of 
all the detector planes with respect to their nominal po- 
sition (global alignment) have been investigated in detail 
to define possible misalignment biases of the track mea- 
surements |17j . The most important internal alignment 
parameters which can be resolved within the T2 hit resolu- 
tion are the shifts of the planes in the X and Y directions. 
Two different methods (iterative and MILLIPEDE [T5] ) 
were used to correct for such displacements. Both gave 
consistent results, with an uncertainty on the transverse 
position of the plane of about 30 ftm. The relative align- 
ment between the two neighbouring quarters was obtained 
using tracks reconstructed in the overlap region. 

The global alignment of the detector is of main impor- 
tance for the present analysis. This was achieved by ex- 
ploiting the symmetric distribution of the track parame- 
ters and the position of the "shadow" of the beam pipe, a 
circular shaped zone of the T2 planes characterised by a 
very low hit rate due to interactions of primary particles 
in the \rj\ = 5.53 beam pipe cone in front of T2. 

The combination of these methods gave, for each quar- 
ter, the X-Y shift with respect to the nominal position 
with a precision of ~ 1 mm and the tilts in the XZ and 
YZ planes with a precision of ~ 0.4 mrad. 

After the local and global alignment parameters had 
been measured with the data, the corresponding misalign- 
ments were introduced into the GEANT4 simulation and 
the same algorithm for the correction of the hit positions 
was applied to the reconstruction of both simulation and 
data. In this way one can also take into account, in the 
simulation, the non-uniform effect that misalignment has 
on the reconstructed hit position, which depends on its 
X-Y coordinate in the GEM plane. 



Data and MC sample. — The sample used for the 
present analysis consists of 150,000 pp collisions at y/s = 7 
TeV recorded in May 2011 during a low pile- up run at 
standard optics. Each beam was composed of six bunches 
with an average luminosity per colliding bunch pair of 
about 8 x 10 27 cm -2 s" 1 corresponding to an inelastic pile- 
up probability of ~ 3%. The rate of beam gas interactions 
for such beam conditions is expected to be negligible. The 
trigger required at least one trigger-road, defined as more 
than 3 "superpads" (3 radial and 5 azimuthal neighbour- 
ing pads) fired in the same r-</> sector of different planes of 
the same T2 quarter. This condition is satisfied if at least 
one charged particle traverses the T2 detector. The ob- 
served trigger rate was about 3 kHz. The fraction of dead 
or noisy channels in this data sample has been measured 
to be 6% for the pads and 9.5% for the strips. 

With the requirement of at least one reconstructed track 
in the T2 detector, the visible cross section seen by T2 
has been estimated to be about 95% of the total inelastic 
cross section. This is based on the comparison of the di- 
rect measurement of the T2 visible inelastic cross section 
(to be published) to the TOTEM inelastic cross section 
measurement deduced from the difference between the to- 
tal and elastic cross sections [20]. The fractions of the 
total inelastic cross section visible to T2 obtained from 
Pythia 6.42 [21] and Pythia 8.108 [22] are in agreement 
with the above estimate. The T2-triggered sample con- 
tains more than 99% of all non-diffractive events and all 
single and double diffractive events having at least one 
diffractive mass larger than ~ 3.4 GeV/c 2 [IB] . 

The transverse momentum (Pt) acceptance for single 
charged particles going into T2 is limited by the magnetic 
field and multiple scattering effects. Simulation studies 
have shown that the charged particle tracks are recon- 
structed, within the analysis cuts utilised in this work, 
with a good efficiency for Pt > 40 MeV/c, defining ef- 
fectively the minimum Pt acceptance. The fraction of 
charged particles with Pt < 40 MeV/c produced in the 
T2 acceptance is predicted to be very small (~ 1%). 

Analysis procedure. — Our pseudorapidity density 
measurement refers to charged particles with a lifetime 
longer than 0.3 x 10 -10 s, and to the charged decay prod- 
ucts of particles with shorter lifetime, which is consistent 
with the ATLAS [23], ALICE [2J] and CMS [25] defini- 
tion of a primary charged particle. With this definition, 
decay products of the K° and A hadrons are considered 
secondary particles, together with all of the charged par- 
ticles generated by interactions with the material in front 
and around T2. The 77- value of a track is defined here 
as the average pseudorapidity of the T2 track hits, calcu- 
lated from the angle that the hit has with respect to the 
beam at the IP. This definition has been adopted on the 
basis of detailed MC simulation studies to find the optimal 
definition of the true 77 of a particle produced at the IP. 
The pseudorapidity density has been measured for each 
quarter independently, allowing an important consistency 
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check among the four analysis results, as each quarter dif- 
fers in its alignment and track reconstruction efficiency. 

Since about 80% of the T2 reconstructed tracks are 
secondaries, it is important to have a procedure for the 
discrimination between them and primary charged parti- 
cles. Based on detailed simulation studies, the most ef- 
fective primary /secondary particle separation is achieved 
using the ZImpact track parameter [17] . This parame- 
ter is proven to be stable against misalignment errors and 
is well described by a double Gaussian function for the 
primary particles and by an exponental function for the 
secondaries. The track ZImpact distribution, with the 
exponential and double Gaussian fit, is shown in fig. [5] 
for data tracks reconstructed in one T2 quarter in the 
5.35 < 77 < 5.4 range. 
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Fig. 2: ZImpact parameter distribution for the data tracks 
reconstructed in one T2 quarter in the 5.35 < 77 < 5.4 range. 
The reported x 2 / n df refers to the global (double Gaussian + 
exponential) fit, performed in the range from -15 m up to 9 
m. The blue-solid curve represents the exponential component 
due to secondaries, while the red-dashed curve is the double 
Gaussian component mainly related to primary tracks. The 
condition Zo • sign (77) < 13.5 m has been required, reducing the 
amount of secondary tracks by about 60%. 

The mean, required to be the same for both Gaussians, 
the standard deviation and the amplitude of the two Gaus- 
sians for primaries as well as the mean and the amplitude 
of the exponential for secondaries have been left free in 
the fit. Since the fit results have been found to be in- 
dependent, the fit, performed on data, has been repeated 



for each 77-bin of the pseudorapidity distribution giving 
standard deviations (amplitudes) of both Gaussians that 
increases (decreases) with 77. The relative abundance of 
secondary particles has been found to be smaller for higher 

The primary tracks were selected using the ZImpact 
and Zo parameters: ZImpact was required to be in the 
range for which the area covered by the double Gaussian 
is 96% of the total, while Zo was required to fulfil the 
condition Zo • sign(?i) < 13.5 m. The fraction of primary 
tracks, among the ones passing the above selection cri- 
teria, was calculated for each 77-bin as a function of the 
ZImpact- value using the double Gaussian and exponential 
fits. This fraction, found to range from about 75% (lower 
1 77 1 bins) to about 90% (higher |?y| bins), allows each data 
track to be weighted by the probability for the track to be 
primary, according to its 77 and ZImpact- value. 

Each track has then been weighted for the primary track 
efficiency according to its 77 and to the pad-cluster multi- 
plicity in the corresponding quarter. This efficiency, eval- 
uated by MC, is defined as the probability to successfully 
reconstruct a GEANT4 generated primary track that tra- 
verses the detector and yields the Zo and ZImpact param- 
eters within the allowed region. The dependence of the 
efficiency from the pad-cluster multiplicity was included 
to make this correction independent from the tuning of 
the MC multiplicity. Once the combined Zo and ZImpact 
requirements were applied, an average primary track effi- 
ciency of ~ 80% was obtained. 

A small contribution to the double Gaussian peak is 
given by the decay products of strange particles and by 
conversion of photons from tt° decays in the material near 
T2. The overall non-primary contribution to the central 
peak, to be subtracted, has been estimated using different 
MC generators and found to range between 6% and 13% 
according to Pythia, Phojet [55] and Sherpa [57]. For im- 
proving the description of the non-primary fraction of the 
tracks, the K° s dN/d?? and 7 dN/dE in the MC were nor- 
malised in order to reproduce, in their acceptance region, 
the measurements by CMS gS] and LHCf [2"9"] . 

MC studies also provided the fraction of primary 
charged particles which do not arrive in the T2 detec- 
tor. The associated correction factor, in average ^1.04, 
has been calculated for each 77 bin in events with at least 
a charged particle in the 5.3 < \rj\ < 6.5 range, by con- 
sidering the number of primary GEANT4 tracks cross- 
ing the detector and the corresponding number of pri- 
mary charged particles generated at the IP with Pt > 40 
MeV/c. 

A bin migration correction accounting for all smearing 
effects on the reconstructed track 77 was also derived. The 
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track distribution was then normalised to the full accep- 
tance in azimuthal angle. The pile-up probability, deter- 
mined from data, was finally subtracted by means of an 
overall correction factor of 0.97. 

Events characterised by a high hit multiplicity, typically 
due to showers generated in interactions with the material, 
were not included in the analysis. These events, where 
track reconstruction capability is limited, constitute about 
11% of the data sample, and have an average pad cluster 
multiplicity per plane larger than 70. The effect of not 
considering these events has been evaluated in a MC study, 
giving an overall correction factor on the measurement of 
about 1.02. 

Eq. (fTJ) was used for the dN c h/d77 determination: 770 is 
the 77- value of the bin centre, S is the sample of tracks with 
r/Q — A77/2 < 77 < 770 + A77/2 satisfying the selection criteria 
above, A77 = 0.05 is the bin width, Ne v is the number of 
events in the data sample, W is the probability for a track 
to be primary, e is the primary track efficiency (where m 
indicates the event pad-cluster multiplicity), Bj is the bin 
migration correction associated with the jth bin, G is the 
correction factor for primary particles not reaching T2, S p 
is the correction factor for the non-primary contribution 
to the double Gaussian peak, $/2tt is the azimuthal accep- 
tance, H is the correction factor taking into account the 
effect of the exclusion of the events with high secondary 
multiplicity and P is the pile-up correction factor. 

Systematic uncertainties. — The systematic uncer- 
tainty associated with the primary track efficiency has 
been evaluated in studies where tracks were reconstructed 
using a set of 5 detector planes (out of the total of 10) in 
a single T2 quarter. The track reconstruction efficiency 
was determined using the other set of 5 detector planes in 
the same quarter. The associated systematic uncertainty, 
estimated with this procedure, was defined as the differ- 
ence between the result obtained using the above data 
driven method, and the MC analysis using the same def- 
initions. This uncertainty, computed as a function of the 
pad-cluster multiplicity and of the track 77, has been found 
to give a relative contribution from 1 to 7 %. 

In order to evaluate the systematics due to the align- 
ment corrections, the global alignment parameters have 
been varied around the optimal values within their reso- 
lution. The data have been again reconstructed and anal- 
ysed for the different misalignment configurations. The 
corresponding variation in the dN c h/d?7 result defined the 
systematics due to the alignment. This uncertainty, which 
is 77 and quarter dependent, has been found to be in the 
range of 3-4%. 

The systematic uncertainty associated with the fraction 
of the non- primary contribution to the central peak, S p , 
has been evaluated by considering the maximum variation 
obtained from several MC generators. It has been found 
to be 77-dependent, ranging from 1% to 3%. 

The systematic uncertainty on the W function, needed 
for the primary to secondary separation, has been esti- 



mated to be in the range of 2-3% considering the uncer- 
tainty on the fitting parameters. 

The primary track efficiency variation due to magnetic 
field effects and to the uncertainty on the energy spectrum 
resulted in an error of about 2%. This contribution has 
been evaluated in the simulation by switching on and off 
the magnetic field and by varying in a reasonable range 
the energy spectrurrQ. 

The variation of the G function from the MC genera- 
tor has been estimated to be around 2%, while the Bj 
functions have been found to have a negligible variation. 
The effect of the track quality criterion requirement, x 2 - 
probability > 1%, has been estimated to be around 1% by 
evaluating the data/MC discrepancy observed with and 
without using this requirement. The effect of the trigger 
bias in our measurement has been evaluated by compar- 
ing the data selected with a pure bunch crossing trigger 
with the sample triggered with T2 and has been found to 
be around 1%. The pile-up probability systematic uncer- 
tainty has been estimated to be 1%. The uncertainty on 
the correction accounting for the exclusion of events with 
high multiplicity of secondary particles has been found 
to be about 1% from the difference between Pythia8 and 
Sherpa MC predictions. 

Table Q] shows the uncertainties of the bin centred at 
770 = 6.025, for one of the T2 quarters. The double dashed 
line separates the quarter dependent contributions (top) 
from the ones in common for all the quarters (bottom). 
The total systematic uncertainty has been computed by 

Table 1: Summary of the relative uncertainties in the bin cen- 
tred at 770 = 6.025 in one of the T2 quarters. The first two 
contributions are quarter dependent. 

770 = 6.025 dN c h/d?7 error summary (one quarter) 



1. Primary track efficiency 4% 

2. Global alignment 3% 

3. Non-primaries in the central peak 2% 

4. Primary to secondary separation 2% 

5. B-field and energy spectrum 2% 

6. Primaries not arriving in T2 2% 

7. Track quality criterion 1% 

8. Trigger bias 1% 

9. Pile-up probability 1% 

10. Events with high secondary multiplicity 1% 

11. Statistical 0.7% 
Total (single quarter measurement) 10% 



1 The CMS experiment reported a discrepancy, that increases with 
r/|. between data and Pythia inelastic event simulation 30 . Pythia 
underestimates the energy flow in the CMS HF calorimeter, mea- 
sured in the pseudorapidity range 3.15 < < 4.9. The effect of this 
discrepancy on the reported analysis has been taken into account us- 
ing a dedicated simulation, where the input energy spectrum of the 
particles has been increased according to the extrapolated discrep- 
ancy expected in the 5.3 < |r7| < 6.5 region. 
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first linearly adding the global alignment and track effi- 
ciency systematics to take into account misalignment ef- 
fects on the primary track efficiency estimation, then this 
result has been added in quadrature to the uncertainty 
contributions from 4) to 10) of tabic Q] and finally the un- 
certainty associated with the non-primary contribution to 
the central peak has been added linearly. To obtain the 
total uncertainty of the single quarter measurement, the 
statistical error is then added in quadrature. 

Results. — The dN c h/d?7 measurements obtained 
for the different T2 quarters are compatible within the 
quarter-dependent systematic uncertainties. For each r\- 
bin, the measurements for the four quarters have been 
combined with a weighted average using only the quarter- 
dependent uncertainties. 

A conservative approach has been adopted for the com- 
bination of the quarter-dependent systematic uncertain- 
ties: an error propagation on the weighted averages has 
been applied, considering the measurements completely 
and positively correlated. The resulting error has then 
been combined with the systematic contributions that are 
common to all quarters and with the statistical one, as in 
the case of the single quarter measurement. 

The pseudorapidity density measurement is shown as 
black squares in fig. [3] where the error bars represent the 
total uncertainty including the statistical error. The pseu- 
dorapidity bins at the edges of the T2 detector, where large 
corrections due to the reduced geometrical acceptance are 
needed, are not reported here, nor are the three bins in the 
5.425 < \r\\ < 5.625 range, where large corrections should 
be applied because of the interaction of the primary parti- 
cles with the beam pipe. Fig. [3] also shows the comparison 
of the data with some MC expectations. Phojet 1.12 (red 
triangles) estimates a ~ 30% (~ 20%) lower dN c h/d?7 than 
measured at \rj\ = 5.3 (6.4). Both Pythia 8.108 with de- 
fault tune (blue circles) and Pythia 6.42 D6T estimate a 
~ 20% (~ 12%) lower dN c h/d?7 than measured at \rj\ = 5.3 
(6.4). Sherpa 1.3.0 with default tune (orange diamonds) is 
in agreement with the data for \rj\ < 5.9, while it estimates 
a higher dN c h/d?7 by ~ 25% at \rj\ = 6.4. Similar trends in 
the comparison of the data with MC predictions have also 
been found by the ATLAS [23J, ALICE [SJ, CMS HO and 
LHCb [3T] experiments in their own r\ acceptance. When 
comparing with cosmic ray MC generators, SYBILL [33] 
predicts a 4-16% lower dN^/dry, while QGSJET01 [33], 
QGSJETII [34] and EPOS [35] predict a 18-30% higher 
dN ch /d7 ? . 

The dN c h/d?7 measurement is also reported in tabled] 
for each 77-bin with the corresponding systematic and sta- 
tistical error. 

Conclusions. — The TOTEM experiment has mea- 
sured the charged particle pseudorapidity distribution in 
pp collisions at y 7 ! = 7 TeV for 5.3 < \ri\ < 6.4 in 
events with at least one reconstructed track in this range. 
This extends the measurements performed by the other 
LHC experiments to this previously unexplored forward rj 
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Fig. 3: Charged particle pseudorapidity density distribution. 
The experimental points (black squares) represent the average 
of the four T2 quarters, with the error bars including both 
statistical and systematic error. Red triangles, blue circles and 
orange diamonds show, respectively, the Phojet, Pythia8 and 
Sherpa predictions for charged particles with Pt > 40 MeV/c 
in events where at least one charged particle is generated in 
the 5.3 < \r)\ < 6.5 range. 



Table 2: TOTEM dN c h/d?7 measurement for inelastic pp 
events at y/s =7 TeV. The reported values represent the aver- 
age for the four T2 quarters, with the corresponding system- 
atic (syst) and statistical (stat) error. 770 represents the central 
pseudorapidity value in each 77 bin. 



Vo 


dN ch /d77 


syst 


stat 


5.375 


3.84 


0.37 


0.01 


5.425 


3.64 


0.38 


0.01 


5.625 


3.54 


0.33 


0.01 


5.675 


3.50 


0.32 


0.01 


5.725 


3.40 


0.30 


0.01 


5.775 


3.42 


0.31 


0.01 


5.825 


3.32 


0.29 


0.01 


5.875 


3.27 


0.29 


0.01 


5.925 


3.20 


0.28 


0.01 


5.975 


3.07 


0.27 


0.01 


6.025 


3.04 


0.26 


0.01 


6.075 


2.94 


0.26 


0.01 


6.125 


2.80 


0.25 


0.01 


6.175 


2.74 


0.26 


0.01 


6.225 


2.65 


0.24 


0.01 


6.275 


2.58 


0.23 


0.01 


6.325 


2.53 


0.21 


0.01 


6.375 


2.38 


0.21 


0.01 
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range. The measurement refers to charged particles with 
P T > 40 MeV/c and with a mean lifetime r > 0.3 x 10~ 10 
s, directly produced in pp interactions or in subsequent 
decays of particles having a shorter lifetime. A prelimi- 
nary measurement of the T2 visible inelastic cross section 
confirms that about 95% of the inelastic pp events have 
been considered in the present study. This comprises more 
than 99% of non-diffractive processes and the single and 
double diffractive processes with diffractive masses above 
~ 3.4 GeV/c 2 . The pseudorapidity density has been found 
to decrease with increasing |?7|, from 3.84 ± O.Ol(stat) ± 
0.37(syst) at |ry| = 5.375 to 2.38 ± O.Ol(stat) ± 0.21(syst) 
at \rj\ = 6.375. Several MC generators have been compared 
to data; none of them has been found to fully describe the 
measurement. 
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